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Abstract: The physics of plasmas containing positive and negative ions is discussed 
with special attention to the recently produced pair-ion plasma containing ions of equal 
mass and opposite charge. The effects of the density gradient in the direction perpendic- 
ular to the ambient magnetic field vector, observed in the experiment, are discussed. In 
addition, the possible presence of electrons is discussed in the context of plasma modes 
propagating at an angle with respect to the magnetic field vector. The electrons may 
either be added to the plasma or enter the plasma attached to negative ions and then 
become released from the ions as a result of collisions. It is shown that the electron 
plasma mode in a cold plasma may become a backward mode in the presence of a den- 
sity gradient, and this behavior may be controlled either by the electron number density 
or the mode number in the perpendicular direction. In plasmas with hot electrons an 
instability may develop, driven by the combination of electron collisions and the density 
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gradient, and in the regime when the ions' response is similar to a sound mode, i.e., 
for un-magnetized ions, implying mode frequencies much above the ion gyro-frequency or 
mode wave-lengths shorter than the ion gyro-radius. In the case of a pure pair-ion plasma, 
for lower frequencies and for parameters close to those used in the recent experiments, 
the perturbed ions may feel the effects of the magnetic field. In this case the plasma 
mode also becomes backward, resembling features of an experimentally observed but yet 
unexplained backward mode. 



PACS Numbers: 52.27.Cm; 52.30.Ex; 52.35.Fp 
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I. INTRODUCTION 



Space plasmas only exceptionally include a single ion species. Much more frequently, 
they consist of a mixture of several ion species with different charges and degrees of 
ionization, different masses and concentrations, on a background of usually more mobile 
electrons. In the laboratory, the situation is not much different. A particularly intriguing 
laboratory example of a multi-ion plasma is the case of the so-called pair-ion plasma. The 
term 'pair-ion' implies a plasma containing oppositely charged ions with nearly the same 
mass. In the recent past, there has been a lot of activity on the theory and the experiments 
dealing with such plasmas, initiated by a series of works describing the experimental 
production of pair ions.^~^ The ions in question are in fact C^, produced in simultaneous 
processes of impact ionization and electron attachment, that are further collected by a 
magnetic filtering effect, i.e., by a diffusion in the radial direction (perpendicular to the 
magnetic field lines). In this way the ions are first separated from the electrons, and then 
collected through an annular hole (with 3 cm diameter), due to a subsequent downstream 
along the magnetic lines, in a very narrow and elongated chamber (90 cm in length, and 
with an 8 cm diameter). In Refs. 2-5 several types of modes have been reported, e.g., 
the ion plasma wave (IPW), the ion acoustic wave (lAW), and the so-called intermediate 
frequency (IF) wave. Detailed measurements presented in the most recent Refs. 4 and 5, 
reveal that the I AW actually has two separate branches, accompanied by some additional 
backward propagating mode situated between them. 

The works 1-5 have attracted a lot of interest from researchers. This has resulted 
in numerous theoretical studies^"^^ dealing with various aspects of linear and nonlinear 
waves and instabilities in pair plasmas. The IF wave has been particularly in focus in 
a number of works, because the behavior of this mode does not fit into the standard 
plasma theory. Yet, the lAW also seems to be problematic because it should not appear 
in the perfectly anti-symmetric pair-ion plasma in question here (i.e., with two ion groups 
with equal masses and temperatures, but opposite charges). However, specific conditions 
in the experiments^"^ related to the different processes (electron impact emission and 
attachment) involved in the production of the positive and negative ions, reveal a small 
disparity in the two temperatures, and this should in principle allow for an electrostatic 
acoustic mode to develop. 

The very recent successful production of pair-hydrogen plasmaP is even more attractive 
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because of the small ion mass and the obvious consequences related to this. This indicates 
that the experimental and analytical efforts in this particular field will increase even 
further. Therefore, in what follows we discuss various aspects of the stabihty of waves 
in pair-ion plasmas, using some limits and parameters that are not necessarily always 
directed only to the presently available experimental results. This study includes various 
limits of the plasma density, and the presence of electrons which may naturally appear 
in such a plasma or can be added for the sake of scientific curiosity to check the effects 
introduced by their presence. 

II. MODEL 

In the experiments^"^ due to the specific filtering mentioned above, only pair ions 
are collected in the chamber, and it is believed that electrons are generally absent. The 
arguments supporting this opinion are convincing and, hence, we accept the mentioned 
filtering as a fact. Yet, electrons may be added to the pair-ion plasma in order to check 
the effect of their presence on the stability of the plasma. On the other hand, we note that 
electrons, being attached to the negative ions, still enter the chamber indirectly. Moreover, 
collisions enable charge exchange between the ions, formally described by a process of the 
type + B~ — > A + B, and resulting in a number of neutrals in the system. The 
electrons can also be released from the negative ions in any of the collisions of the type 
B-' + B B + B + e~, B- + A ^ B + A + e~, SindA+ + B~ A+ + B + e~. As a result, 
a mixture of pair-ions (because in the present case A = B = Ceo), electrons, and neutrals 
may be expected. Hence, the momentum equations for ions and electrons, that will be 
used below, will include the friction force between different plasma species. According to 
the given scenario, the presence of electrons should be related to the presence of neutrals. 
Yet, the neutrals are not expected to play a significant role in the system because in 
terms of electron collisions, roughly speaking, an ion is equivalent to 3.4 ■ 10^(300/Te)^ 
neutral atoms. So collisions between charged species are usually dominant even when the 
number density of neutrals n„ by far exceeds the number density of ions. The situation 
is quite opposite with the presence of electrons. As a matter of fact, even a small amount 
of electrons can drastically change the dispersive properties of the modes in multi-ion 
plasmas, and indeed introduce new ones. 

We note also that the spatial variation of the density in both the radial and the axial di- 
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rections has been experimentally detected/ with the radial density profile being Gaussian. 
In view of the very different axial and radial lengths of the chamber, the axial density 
inhomogeneity may be neglected without losing any essential effect. Such equilibrium 
density gradients are a rather common feature of various laboratory plasmas. ^''"^^ 

In the experiments 2-4 the modes are initiated by a cylindric exciter at one end of the 
chamber, and they also appear to propagate almost parallel to the magnetic field vector, 
i.e., they may include a (small) wave number component in the perpendicular direction 
as well, as commented in Ref. 4. 

A. Unmagnetized ions 

We proceed by introducing high frequency electrostatic perturbations propagating 
obliquely to the equilibrium magnetic field vector Bq = BqCz. In the frequency limit 
Qe > ^ ^i, the ions are un-magnetized, the Lorentz force in the ion momentum 
equation can be omitted, and their perpendicular (with respect to the magnetic field 
vector) response to the perturbed electric field is similar to the case of an ordinary ion 
sound wave. In principle, the ions will remain un-magnetized even for frequencies below 
the ion gyro-frequency (that are in fact also experimentally measured) provided that the 
perpendicular wave-length is below the ion gyro-radius. The momentum equation for ions 
of type j is as follows 

dv ■ 

^i^j-gf- = -^jlj^^ - i^TjVrij - mjUjUjVj. (1) 

The electron dynamics is described by the momentum equation 

dv, 



at 



eneV0 — eUeVe X B — KTeVrie — rrieneiyeVe, (2) 



and we also use the appropriate continuity equation for all species 



dt 



+ V ■ {rijVj) = 0. (3) 



'en- 



Here, the electron collision frequency is in principle the sum Ue = fei + ^'e 

From Eqs. ([1]) and ([3]) we get,^° for ion perturbations of the form ~ f{x) exp{—iujt + 
ik ■ r), where /^denotes an arbitrary direction with respect to Bq, 

9 I ■ 19 2 1 

rrij (jj^ + lujUi — /e^fyj 
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where v'^j = KTj/rrij. The electron dynamics in the perpendicular direction is described 
by the following recurrent formula for the velocity 

Ve± = -=r(^z X V±0 + -- — — X 



He ni Tie 'T-e 

-^(| + .^e-v)e.x.-.^-^(| + .- .v).-ex. (5) 
The parallel perturbed electron velocity is given by 

Vezl = 01 + — ■ (6) 

rUeUJl Ui UeO 

Here, ui = Uq + iu^, uq = u — kyVeo, and Veo ~ — (e^ x V _LneQ/neo)v'le/Qe is the electron 
diamagnetic drift. We have used u^, <^ fig in the electron diamagnetic drift and in the last 
term in Eq. and we work in the usual drift limit for electrons, with a small equilibrium 
density gradient, and in a local approximation implying \d/dx\ <^ \ky\. 

III. MODES IN COLD PLASMA 

In order to check for the most basic modes, here we omit collisions and thermal effects, 
and discuss first the case of a two-component electron-ion plasma. The simplified ion equa- 
tion (i.e. without thermal terms, Hj) is used, while the electron density perturbation be- 
comes Uei = [ueoky/ (QeBo) — CHeokl / {nieUj'^) — kyn'^^^/ (c<ji?o)]0i. From the quasi-neutrality 
condition we obtain a modified lower hybrid mode, viz. 

0_ (7) 

kyHeO k^ \ ilekij 

Using the Poisson equation instead of the quasi-neutrality, we obtain 

^'fl + ?i51-T#^--^p-^ = 0. (8) 
\^ Q^k^ J k^Qerieo ^ k^ 

Comparing this to Eq. ([7]), it is seen that the quasi-neutrality in the perturbed state 
is equivalent to assuming Up^k"^/ {Qlk"^) ^ 1. For perpendicular propagation without 
a density gradient, Eqs. (I7l) and (IE]) yield the lower- hybrid and the ion plasma waves 
= Q^Q. anfi ^ ujpi/{l + cUpg/fig), respectively. 

In view of the difference in mass, the last three terms in Eq. ([8]) can not always be of the 
same order. An interesting case for which the mode behavior is determined by the density 
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gradient term, yielding a backward solution du /dk^ < 0, is obtained by the following set 
of parameters. We take Bq = 0.1 T, L„ = rieo/n'^Q, where neo{r) = NQexp{—r^/a^), so 
that locally = a^/(2r). If a = tq, or a = 2ro, where ro = 4 cm is the plasma radius, 
at r = ro this gives for the density 0.37 and 0.78 times its value at the axis, respectively. 
We choose L„ = 8 mm, and take three values for the plasma number density, omitting 
the index e for a quasi-neutral equilibrium density: Uq = 10^^ m~'^, Uq = 5 ■ 10^^ m~^, 
and uq = 10^^ m~^. The results are presented in Fig 1. Here, the wave frequency has a 
value much above the ion gyro frequency, but below the electron gyro frequency, and the 
parallel wave length Xz is in the range 10 — 100 cm. In the poloidal y (or 6) direction we 
have fixed one wave length only, i.e., m = 1. The decreasing (full) line with (backward 
mode) represents the solution of Eq. ([H]) for the case of a dominant density gradient term. 
For a lower plasma density (dotted line) we have a direct mode du/dk^ > 0. In a cylindric 
geometry, the given solutions describe twisted waves, i.e., those traveling both along the 
magnetic field lines and in the ^-direction. The backward mode has a behavior similar to 
the observed IF wave, yet in the much higher frequency domain. 

The same effect, i.e., the transition from a backward to a direct mode, can be achieved 
and controlled by changing the poloidal wave length Aj,. This can easily be checked by 
setting uq = 10^^ m~^, and taking m = 10, which yields a direct mode. 

It is worth mentioning that this transition from forward to backward modes by increas- 
ing the pressure, has in fact been observed experimentally in rare gases as was reported 
long ago in Ref. 21, but without any explanation about the source of that behavior. 

We stress also that in the presence of an additional, negatively charged ion species of 
the same mass, the mode behavior does not change. This is because the case presented 
here is an interplay between the two dominant terms, i.e., the density gradient term and 
the last term in Eq. ([8]). 

IV. GROWING MODES IN COLLISIONAL INHOMOGENEOUS PLASMA 



Using Eqs. ([5]) & ([6]), from the electron continuity equation we obtain: 

riel ^ ^* - vleKi^ - aUJiUJak'^/{klQl)]/uJi 601 

rieo uj -vlJvl[l- auJiUJaky{klVtl)]/uji kT^' 

Here, Ua = UQa + iUe, « = Vll + '^e/^e)) ^* — ~kyKTen'^Q/ (cBQUeo). Note that Uq appears 
in both uJa and uji from the left-hand side of the electron momentum equation, i.e., as the 
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finite electron mass effect. 



A. Massless hot electrons 



22 



In the limit |a;o| < i^e, from Eq. ([9]) we have 

The second term in Dg is small compared to the first one provided that kyP^/{k'^fll) -C 1, 
implying a dominant effect of collisions on the parallel electron dynamics. In this case, 
Eq. (fTOl) coincides with the corresponding expression from Ref. 23. The same limit holds 
when the electron mass corrections are retained. ^^'^"^ 

Using Eq. (fTOl) instead of Eq. ([9]), with the quasi-neutrality in the equilibrium and the 
Poisson equation in the perturbed state, in the case of the pair-ion electron plasma with 
Qa = —qb = e, we obtain the dispersion equation 



Here, u^^ = e^rijo/ {sQUij). Below we discuss several particular cases. 



a) Without electrons and for a collision-less, pure pair- ion plasma with = m;,, and 
for modes propagating along the magnetic field vector, Eq. ffTTl) yields 



^2 _ ^2 _|_ ^li'^Ta + "^Tfc) _j_ 



1/2 

(12) 



In the limit of a small kz (and/or a small difference between the two temperatures) 
this describes the ion plasma wave uj"^ ~ 2^01^ + kliy"^^ + v\^)l1, and the ion sound 
mode up' ~ k^v^^ + vlf^)/2 when 7^ Tf,. It is tempting to set = here, 
however, strictly speaking, the sound mode loses its electrostatic nature in this case 
and it becomes an ordinary gas acoustic mode (or thermal mode) involving only the 
pressure perturbations. 

Specific conditions in the experiments,^"^ related to different processes involved in 
the production of positive and negative ions (electron impact emission and attach- 
ment), reveal a small disparity in the two temperatures Ta 7^ Tb. Hence, solving 
Eq. (fT2l) with the parameters riio = 10^^ m^^, Ta = 0.5 eV, Tb = 0.3 eV, and for kz 
in the range 6 — 125 m~^, yields the ion plasma wave with frequency 220 kHz, and 
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the ion sound wave in the range 1.5—29 kHz. The sound mode in the experiments^"^ 
appears divided into two branches. The values for the ion sound which follow from 
Eq. (fT2l) are in good agreement with the lower experimental branch. However, for 
explaining the upper experimental branch the presence of electrons may be needed. 
In view of a large difference in mass, the electrons are normally described by the 
Boltzmann distribution. Yet, adding such a contribution of electrons only yields 
corrections that are of importance at spatial scales of the order of the electron De- 
bye scale Xde = VTe/^pe- Therefore, extra effects are needed, like those that follow 
from Eq. (I9l), and from the last term in (ITTll . 

b) In an electron-ion (nj, = 0) quasi-neutral plasma (A;^A^g ^ 1) with cold ions, from 
Eq. ffTTj) we have 



describing a high-frequency (as assumed u > VLi) modified ion acoustic mode that 
may be unstable due to the density gradient and electron collisions. The numerical 
solution for the growth rate of the strongly growing mode fll3p is presented in Fig. 2, 
for a relatively high (necessary to have all conditions satisfied) plasma number den- 
sity rieQ = 10^^ m"'^, the magnetic field Bq = 0.2 T, for cold ions, and = 5 eV, 
L„ = 2 cm, and for the three values ~ 1, 2, 5 mm. For these parameters and 
the given A^,, we have u^, = 7.8 ■ 10^, 4.7 ■ 10^, 1.6 ■ 10^ Hz, and the corresponding 
modified sound wave frequency ~ 5 ■ 10^, 3 ■ 10^, 1 ■ 10^ Hz, respectively. The 
mode is growing for uj^ > uo. Here mj = 720 mp, hence f2j <^ and the ions are 
un-magnetized. In the same time, fc^A^^ ^ 1 for all three cases discussed, so that 
the quasi- neutrality is well satisfied, while tuo/z^e ^ 1 so that the electron inertia is 
negligible. The FEM (finite electron mass) line is described further in the text. 

c) In the case of pair-ion electron plasma with the Boltzmann electron response (that 
implies ^ oj and u/k <^ Vtc), the last term in Eq. (ITTl) is to be replaced with 
cUpg/ (A;^f|,g). Setting Ta = Tf, = T and nia = nib = m yields only one pair of solutions 



00^ = ecl[T/T, + {nao + nto)/n,o] for PX^ « 1, and = (cu^^ + ^1^)^^X1/(1 + 
k'^Xl) + k'^v^ for fc^Ag > 1, where = nTe/m, = kT/itl. 



Otherwise, from Eq. (ITTl) for cold ions we have the third order dispersion equation 




(13) 




(14) 
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Here, s = l + l/(fc^A^g). Unstable solutions may be detected by using the generalized 
Hurwitz method for polynomials with complex coefficients. According to this, 
for a polynomial of the degree m, of the form (ao + ibo)x^ + (ai + ibi)x 



m—l 



+ . . . + 



{am + ibm) = 0, one makes the sequence of m + 1 numbers cq = ao, ci 
. . ., where r goes to m, and where 



ai. 



■1) 



r(r-l)/2 



ai 


ao 











-b2 


-61 


ai 


ao 





as 






&i 


ai 




-6s 


as 


02 


-62 ■ 


ag 


a4 




bs 


as 


0-2r-l 


0'2r~2 


&2r-2 


&2r-3 


02r-3 ■ 



The number of roots with positive real parts equals the number of sign changes in the 
sequence cj. A sufficient instability condition is that any of the Cr has a negative sign. 



Applying this to Eq. ([HD, we find: Cq = 1, ci = u^/{k'^\lJ, C2 = -u'^^uj'^ / {k"^ X^J'^ , 
Cs = —D'^ulup^/^k^Xj^J. This indicates that, for cu* > 0, Eq. (IT^ has one mode 
propagating in the positive direction, and two solutions propagating in the opposite 
direction, and there exists at least one growing solution. Another possibility is 
uj^ < 0. According to experimental conditions n'^Q < 0. Therefore in this case we 
have ky < 0, there are two sign changes in the series cj, hence two positive solutions 
for LJ, and we still have at least one growing solution because ci, C2 < 0. 

d) The case a) does not change considerably in the presence of boltzmannian electrons 
in the collision-less plasma and for the same parallel propagation when we have 



2s 2 ^ ™ 



^Tb) 



1 

±- 

2 



^Tb) 



"^{^pa ~ ^pb) , , 2/ 2 

h kJv, 



1 ^ V2 



Ta 



^Tb) 



This is just a modified form of Eq. (fT2!) . 

B. Full dispersion equation 

Keeping the full form of Eq. ([9]), with the help of the Poisson equation, the dispersion 
equation becomes 



k'^v. 



2 

Ta 



Tb 
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Eq. (|T5l) will be solved numerically for some sets of parameter values. 

We first check the effects of the finite electron mass (FEM) , and solve Eq. ( fT5|) for the 
electron-ion case (nfeo = 0) with the same parameters as in Fig. 2, for cold ions, and for 
the case ky = 3770 m^^ (dashed line). The sound wave frequency appears to be slightly 
larger, i.e., in the range 3.1 — 3.2 MHz, while the growth rate is reduced, given by the 
lower dashed line in Fig. 2. 

For the pair-ion electron plasma, Eq. ( ITSl) is solved for collision-less ions in terms of 
nbo/nao, and the result for the acoustic mode is presented in Fig. 3 for the same parameters 
as in Fig. 2 and for fixed values of kz = 18.8 m~^ and ky = 3770 m^^. The driving source 
for the instability is the density gradient in combination with electron collisions. Hence, 
for a reduced electron number density the instability vanishes; in the present case this 
happens at nbo/uao — 0.4. The other pair of solutions, that is in the range of the ion 
plasma frequency and the lower hybrid frequency, appears strongly damped, with Ur 
being of the same order as \uji\ and it is consequently of no interest. 

For much lower density, e.g., setting n^o = 10^^ m~^ and varying riho/nao, we obtain 
two direct modes with frequencies close to 10^ Hz and 2 ■ 10^ Hz, propagating almost 
without change, the former very weakly growing and the latter being almost undamped. 

In the case of different temperatures of two pair-ion species (like in the reported cases 
in Refs. 2-5), there appears an additional acoustic mode. This is checked by setting 
Ta = 0.5 eV, Tfe = 0.3 eV for Uao = 10^^ m~^. The other parameters are the same as 
in Fig. 3. In total, there are three pairs of solutions: one high frequency ~ 10'' Hz but 
highly damped mode mentioned above, the slightly modified acoustic mode from Fig. 3, 
and the third, weakly growing acoustic mode presented in Fig. 4. 

We stress that the negative solutions (that are not presented in the figures above) are 
in fact not counterpart solutions in the strict sense because of the asymmetry caused by 
the density gradient. 

V. WEAKLY MAGNETIZED IONS 

In the experiments^"^ the ions in the equilibrium are magnetized. However, for per- 
turbations with frequencies above the ion gyro-frequency, the ions are un- magnetized, yet. 
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for not so high frequency oscillations they can feel the effects of the magnetic field. Hence, 
in Eq. ([1]) we add the Lorentz force. The equilibrium diamagnetic drift velocity for the 
two ion species is given by Vjo = {cz x Vnjo/njo)KTj/{qjBo). For the given experimental 
parameters (T^ = 0.5 eV, Tb = 0.3 eV) and assuming L„ of the order of 1 cm, the two 
velocities are of the order of 10^ m/s. For perpendicular wavelengths of the order of 1 cm, 
the Doppler shift in the frequency, kyVjo, becomes of the order 10^ Hz, and can therefore 
introduce visible effects in the mode behavior. Assuming a negative density gradient, we 
have a negative Vaoe^ and Vbo = —VaoTb/Ta- The velocity components for the ion species 
a become: 

Vaxl 



1 



^2 



-Df 



Vayl 



30 



kyUJaO^^ + via 
-DO 



kyQ — UJaO 













nao_ 




\ rial 




) nao 



Vazl 



+ 



kzVia rial 



^aQ riaO 

Here, Uao = uj — kyVao- In the corresponding equations for the ion species 6, the term Vt 
is to be taken with the minus sign. 

The continuity equation Uji/rijo = {kyVjyi + k^Vj^i + iVjxi\nQ/nQ\)/ujjQ for the two 
species, and the Poisson equation, yield the dispersion equation for the modes in an 
inhomogeneous pair-ion plasma: 

K , k^ 



UJ. 



aO 



UJ 



aO 



n 



oO 



^ao(^aO - 



UJt 



k'^v^ 
-L 2 — 



UJ. 



k'^v^ 

S Ta 



aO 



Q 1 

UJaO k,. 



n 



Ml 



riao 



+ 



n 



riao 



V, 



UJ. 



aO 



Ta ^ f 

\ ^aO 



n 



an. 



nao 



+ 



^^^^^ 



UJ. 



+ 



bO 



n 



aO 



riao 



2 



uJbojuJbo — ^ 



k^v"^ 

1 "'z '^Tb 
i ri — 



UJ. 



bO 



'^y^Tb 



n 



nao 



n 



MI 



nao 



_J!Tb 



UJbQ ^ 



n 



riao 



(16) 



In the absence of electrons, the density gradients and number densities are the same, 
and ujp = e^nQl{eQm). Various limits of Eq. (fT6|) can easily be discussed, yielding the 
previously discussed modes. For example, in the absence of a density gradient and for 
strictly perpendicularly propagating oscillations in a quasi-neutral plasma, from Eq. (fT6l) 
we have just the ion cyclotron mode uj"^ = VL^ + k'^i^ra + v\b)l'^- For strictly parallel 
propagation in a quasi-neutral plasma we have the acoustic mode uP' = k'^iv'^a + "^Tfe)) 
otherwise it yields Eq. f|T2|) . 
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A particularly interesting limit is seen also in the case of a cold plasma and for long 
wavelengths, when the remaining density gradient terms in the numerators in Eq. f|T6l) 
exactly cancel each other, and it yields electrostatic convective cells lo^ = Vt^k1/{ky + k1). 
This result is completely different compared to the standard electron-ion plasma which in 
such a geometry has the drift wave driven by the density gradient. 

Eq. flTB]) is solved numerically for parameters close to the experimental values. We 
take Bq = 0.2 T, = 0.5 eV, = 0.3 eV, but with a lower number density n^o = UbQ = 
10^^ m~^, and we take L„ = 0.8 cm, ky/{2'K) = 0.05 cm~^. The result for modes with 
positive frequency is presented in Fig. 5. For values of kz in the range 0.01 — 1.15 (in the 
given units), the ion plasma mode is in fact a backward one (frequency decreases from 
14.05 at kz = 0.01 to 12.7 at k^ = 0.15), and these values are only slightly below the 
experimental values that correspond to the IF mode. However, we can obtain exactly the 
same frequency as for the observed IF mode by simply slightly increasing the density to 
the value n^o = n^Q = 3 ■ 10^^ m~^, but then the decrease with k^ in this case becomes less 
pronounced, though it still remains there. We suspect that the unexplained experimental 
IF mode may in fact be the ion plasma mode presented here, while the experimental 
upper mode (which in Refs. has been interpreted as the ion plasma mode) may be 
an eventual contribution of electrons, described by our Eqs. and flTU]) . However, these 
statements are speculative as we have not been able to obtain an additional electron mode 
in the given experimental range, i.e., at frequencies about 30 and higher (in the units given 
in Fig. 5). 

Due to the density gradient there is a coupling between modes b and c, in the range 
kz ~ 0.032 — 0.07, with the maximum growth rate Ui of the order uoi/uoj. ^ 0.2. The 
corresponding growth rate, multiplied by 10, is presented by dashed line in Fig. 5. The 
density gradient modifies also the acoustic mode d so that for small values kz its frequency 
becomes negative. This all can easily be demonstrated by switching of the density gradient 
terms. 

We note that the line h has a range of kz for which it is backward duj/dkz < 0, and 
the same holds for the line c. The frequency of both modes, in the range where they have 
backward properties, is in the domain of the low frequency backward mode observed in 
the experiment (around 4 in the given units). Though, the agreement with experimental 
solutions is far from perfect, and the splitting of the lines h and c we obtain for kz > 0.07 
has indeed not been observed. 
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VI. CONCLUSIONS 



The pair-ion plasma containing heavy fuUerene ions (with equal positive and negative 
charge), that has been created in the laboratory, has attracted a lot of interest in the 
recent past. The present study deals with some additional properties of such plasmas, 
e.g., those that result from the eventual presence of electrons, from plasma inhomogeneity 
combined with electron collisions, and from much higher densities. 

The arguments given in Sec. 11 show that the presence of electrons can not be excluded, 
or that they can be added to the system in order to check the consequences of their 
presence. In a recent study^ some properties of waves in pair ion plasma have been 
discussed, showing that in fact they can be used as an indicator of the presence of electrons. 
Some indications in this direction can also be found in the most recent experimental 
results,^'^ where the measured sound branch in fact splits into two sub-branches. In the 
present work, the effects of the density gradient and electron collisions in a relatively 
high density pair-ion plasma have been explored showing the presence of a backward- 
propagating high-frequency mode with features similar to the observed IF mode that 
remained unexplained so far. However, our first backward mode, determined by the 
density gradient (see Fig. 1), appears in a much higher frequency range and it is therefore 
unable to explain the observed IF mode. The case of a pure pair-ion plasma, discussed 
in Sec. V, reveals that the ion plasma wave behaves as a backward mode in a large range 
of the parallel wave number, and perhaps may be interpreted as the IF mode observed 
in the experiment. Yet, the agreement with experimental modes is not good enough to 
draw a definite conclusion. Nevertheless, in our view, the plasma features discussed here 
represent a step forward in the theory of such a newly created plasma. In the same time 
we believe that the presented results should make a good basis for future investigations 
directed towards the exact explanation of the IF mode. 
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Figure captions: 



1. The positive solution of Eq. ([8]) for three number densities in electron- ion plasma, 

showing the transition from direct to backward mode. The parameters are given in 
the text. 

2. The growth rate of the sound mode from Eq. (fT3l) for three values of the perpendicular 

wave lengths. FEM describes the difference due to electron inertia. 

3. The positive solution of Eq. f[T^ for a modified ion acoustic mode as a function of 

the negative-ion number density, and for equal temperatures of the two ion species. 

4. Additional acoustic mode from Eq. f|T5l) in the case of different temperatures of 

positive and negative ions, for Ta = 0.5 eV, and Tb = 0.3 eV. 

5. The positive numerical solutions of Eq. ( |T6l) for a pure pair- ion plasma. The growth 

rate of the coupled modes (multiplied by 10) is given by the dashed line. 
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